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Abstract: We analyze the implications of several recent cosmological and experi- 
mental measurements for the mass spectra of the Constrained MSSM (CMSSM). We 
compute the relic abundance of the neutralino and compare the new cosmologically 
expected and excluded mass ranges with those ruled out by the final LEP bounds on 
the lightest chargino and Higgs masses, with those excluded by current experimental 
values of BR(i? — > X s j), and with those favored by the recent measurement of the 
anomalous magnetic moment of the muon. We find that for tan /5 < 45 there remains 
relatively little room for the mass spectra to be consistent with the interplay of the 
several constraints. On the other hand, at larger values of tan (5 the decreasing mass 
of the pseudoscalar Higgs gives rise to a wide resonance in the neutralino WIMP 
pair-annihilation, whose position depends on the ratio of top and bottom quark 
masses. As a consequence, the cosmologically expected regions consistent with other 
constraints often grow significantly and generally shift towards superpartner masses 
in the TeV range. 
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1. Introduction 

Cosmology provides an important restriction on otherwise allowed supersymmetric 
mass spectra. This in particular is the case with the relic density of the lightest 
supersymmetric particle (LSP) which, in the presence of i?-parity, is stable. A natural 
candidate for the LSP is the lightest neutralino \. Its relic abundance Q x h 2 has 
been a subject of a large volume of papers |T|, starting from Refs. [TJ, [| up to the 
recent comprehensive studies {|, ^, || [Tfl. It is well known that VL x h 2 can vary over 
several orders of magnitude but that it is also often consistent with the abundance 
of non-baryonic cold dark matter (CDM) in the Universe whose determinations have 
been improving steadily. Recent reviews give more narrow ranges for the components 
(matter and baryonic f2b) of the Universe and the Hubble parameter h than in the 
past. For example, for the matter component in Ref. || one finds Qm — 0.33 ±0.035 
which is consistent with 0.4 ±0.1 of Ref. § and 0.35 ± 0.05 of Ref. pf. The 



baryonic component is now Q^h 2 ~ 0.02 ^, |H]] while the Hubble parameter is 
h = 0.72 ±0.07 H (0.65 ±0.08 0). Based on this, and assuming that the neutralino 
LSP makes up the dominant component of the CDM in the Universe, we now select 
the range of the neutralino relic abundance 

0-l<^ 2 <0.2 (1.1) 

as conservatively matching the recent observations. Furthermore, we put an upper 
bound 

n x h 2 <0.3 (1.2) 
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based solely on the constraints on the age of the Universe. These ranges will also 
allow a comparison with the favored values of Q x h 2 in the range 0.1 <Q x h 2 < 0.3 
which have often been used in the recent literature. 

The latest determinations of ^cdm from the measurements of the CMBR, based 
on assuming reasonable ranges for other cosmo logical parameters (priors) , like I^tot = 



1, etc, give ^cdm^ — 0.14 ± 0.04 111, O, or, assuming more priors, even much 



smaller errors: f2cDM^ 2 — 0.13 ± 0.01 [|T^]. We do not feel yet ready to accept 
these narrow ranges as robust enough for our analysis. In particular, we note that, 
assuming h = 0.72 ± 0.08, the DASI analysis [0] finds fi M = 0.40 ± 0.15 which 
implies 0.1 < ^m^ 2 < 0.35. Likewise, the most recent ROSAT measurement gives 
= 0.35 ± 0.12 [ 13 1 which again allows for larger values of flyih 2 . Neverthe- 



less, below we will discuss the impact on our results of assuming, after subtract- 
ing the baryonic component, that f2 x /i 2 is in the range 0.10 < Q x h 2 < 0.12 and 
0.10 < tt x h 2 < 0.15. 

Cosmological constraints ( |l.l| )-( |i~2| ) have a particularly strong effect in the frame- 
work of the Constrained MSSM (CMSSM) 0, |f|. In the CMSSM, in addition to the 
requirement of a common gaugino mass m^j at the unification scale Mqut, which is 
usually made in the more generic Minimal Supersymmetric Standard Model (MSSM), 
one further assumes that the soft masses of all scalars (sfermion and Higgs) are equal 
to m at M GU t, and analogously that the trilinear soft terms unify at M GU t at some 
common value Aq. These parameters are run using their respective Renormalization 
Group Equations (RGEs) from Mgut to some appropriately chosen low-energy scale 
Qo where the Higgs potential (including full one-loop corrections) is minimized while 
keeping the usual ratio tan/3 of the Higgs VEVs fixed. The Higgs/higgsino mass 
parameter [i and the bilinear soft mass term Bfi are next computed from the con- 
ditions of radiative electroweak symmetry breaking (EWSB), and so are the Higgs 
and superpartner masses. The CMSSM thus has a priori only the usual 

tan/?, mi/2, m , A , sgn(/x) (1.3) 

as input parameters. (A comprehensive set of formulae that we will refer to can 
be found for example in Ref. [|HJ.) However, in the case of large m 1 / 2 , m ^ 1 TeV 
and/or large tan/5 ~ Oim t jm^) some resulting masses will in general be highly 
sensitive to the assumed physical masses of the top and the bottom (as well as the 
tau) M but they will also strongly depend on the correct choice of the scale Q . This 
in particular will affect the impact of the cosmological constraints fll.l|) -( |i~2"|) as we 
will discuss below. 



In the CMSSM, the LSP neutralino is often a nearly pure bino |L7|, [18], HH because 
the requirement of radiative EWSB typically gives 3> M\ where M\ is the soft 
mass of the bino. This often (albeit not always! |]15|) allows one to impose strong 
constraints from £l x h 2 < 0(1) on m^j and m (and therefore also on heaviest Higgs 
and superpartner masses) in the ballpark of 1 TeV. This was originally shown in 
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Refs. [ 18] , |i~5| and later confirmed by many subsequent studies starting from |L9| up 
to the most recent analyses [j|, |^, |7j . 



2. Experimental Constraints 

In the case of the CMSSM, the most important experimental constraints from LEP 
are those on the masses of the lightest chargino xt an d Higgs boson h. For the 
first one we adopt the bound m ± > 104 GeV since the actual limits from the LEP 
experiments are very close to this value. The lightest Higgs mass has been a subject 
of much debate on both the experimental |21| and the theory side. The termination 



of LEP has left the first one unresolved for at least a few years. On the theory front, 
due to large radiative corrections, the precise value of still remains somewhat 



dependent on the procedure of computing it [|^, ^3], |24|, [25fl . In our analysis we 
will conservatively assume rrih > 113 GeV but will keep in mind that the theoretical 
uncertainty in rrih in the CMSSM is probably of the order of 2-3 GeV. 

Non-accelerator experimental results are also of much importance. First, there 
has been much recent activity in determining BK(B — > X s j). A recent combined 
experimental result [H 1 BR(B X s -f) = (3.11 ± 0.39) x 10" 4 allows for some, but 
not much, room for contributions from SUSY when one compares it with the updated 
prediction for the Standard Model (SM) BR(B -> X s7 ) = (3.73 ± 0.30) x 10~ 4 (SJ. 
Second, at large tan/3 next-to-leading order supersymmetric corrections to b — > 37 



become important |28|, |29], [31], |32], [33] . In our analysis we adopt the full expressions 
for the dominant terms derived in Ref. [Q. We also include the 6-quark mass effect 
on the SM value which was subsequently pointed out in Ref. |ft3| . We add the two 



la errors (the experimental and SM) in quadrature and further add linearly 0.2 to 
accommodate the theoretical uncertainty in SUSY contributions which is roughly 
5% of the SM value for branching ratio [j35|l . Altogether we conservatively allow our 



results to be in the range BR(£> — > X s ^) = (3.11±0.69) x 10~ 4 for SM plus two-Higgs 
doublets plus superpartner contribution. The excluded regions of SUSY masses will 
not however be extremely sensitive to the choice of these error bars. 

Lastly, much excitement has recently been caused by the first measurement by 
the Brookhaven experiment E821 of the anomalous magnetic moment of the muon 
<V = {g^ — 2)/2 P6| . Taken at face value, the result implies a 2.6a discrepancy 
between the experimental value and the SM prediction a^ xpt — a^ M = (43 ± 16) x 
10~ 10 . There is much ongoing debate about improving the understanding of the 
precise contribution from the SM, as has been reported for example in Ref. [[37] . 
In particular, there has been a tendency of moving the SM prediction towards the 
range of experimental values ||38|| . In our analysis we will use the published results 



1 The most recent update [Q, which incorporates the new CLEO result ^7j, gives BR(_B — * 
X s j) = (3.23 ± 0.42) x 10~ 4 . As we will coment again later, adopting this new data would exclude 
somewhat larger regions of the CMSSM parameter space. 
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to allow SUSY contributions in the ranges 27 x 10~ 10 < Ac^ USY < 59 x 10~ 10 (la) 
and 11 x 1(T 10 < Aa^ USY < 75 x 1(T 10 (2a) but will comment on the effect of 
varying these numbers later. Here we only note that we consider the upper limit on 
Aa^ USY (which implies a lower limit on and m ) as rather robust. In contrast, 
the lower limit on Aa^ USY (and the resulting upper limit on mi/2 and mo) should be 
approached with much caution. We will comment on this further when we present 
our results. 



3. Procedure 



We calculate superpartner and Higgs mass spectra using the package ISASUGRA (v. 7. 51) 
but make some important modifications which will be described below. We refer the 
reader to Ref. 



39 



for a more detailed description of the code. Here we will only 
highlight the main points of the procedure. The overall strategy is to first find an 
approximate spectrum of the Higgs and SUSY masses and then iterate the procedure 
of running the RGEs between Mqut and the low scale until a satisfactory consistency 
is achieved. In the initial step one runs the two-loop RGEs for the gauge couplings 
up and identifies Mqut as the point where ai and cti are equal. Exact unification of 
a s with the other two SM gauge couplings is not assumed as it would predict a value 
far too large compared with the experimental range a s (mz) = 0.1185 ± 0.002 [|D| 
which we take as input. 

The top, bottom and tau masses must also be treated with care since their 
assumed values often have an important effect on the running of the RGEs, especially 
at large tan (3. The pole mass of the top m 



pole _ ]^y4 3_)_5 2 j s initially converted 
to the running mass m t (m t )^ in the MS scheme using a two-loop QCD correction 
with a s running computed at three loops, and is then identified with mt(mt)sM m 
the DR scheme. In subsequent iterations, once SUSY masses have been computed, 
also one-loop SUSY corrections from squark and gluino contributions are included 
to give m t (mt)susY- 

As regards the mass of the bottom, the Particle Data Book |40| quotes the 
range 4.0-4.4 GeV for the running mass in the MS scheme m^m^)!^. We adopt 
the similar range m^mi,)^ = 4.25 ± 0.15 [(41]] which has been used in other recent 
analyses [^|, |5J. In ISASUGRA m b (m b )^ is initially run up to m z using a one-loop 
QCD+QED formula, with a s running computed at three loops and a em at one loop, 
to obtain mb(mz)sM which is assumed to be the same as in the DR-scheme. In 
subsequent steps, once the SUSY masses have been computed, one includes SUSY 



corrections |2g, [42| from squark/gluino and squark/chargino loops which become 
important at large tan f3 



m b(mz)svsv 



m b(m z )sM 



SUSY' 



(3.1) 



4 



where m b (mz)sM * s computed in the initial step and (Am b /m b ) SVSY stands for a 
SUSY mass correction. (We note that in ISASUGRA it is the pole mass m b ole that 
is used as an independent parameter and is actually 'hardwired' at m b ole = 4.9 GeV 
[and in some routines at ml° le = 5.0 GeV]. We have modified the package to allow 
it to treat ml° le as an external parameter which is then adjusted to match the 
desired value for m^mf,)^. The conversion between the two is done using a one- 
loop expression.) 

Finally, the pole mass of the r-lepton is now well measured m^° le = 1777.0 ± 



0.3 MeV pO[ . It is converted to the DR scheme in SUSY using an analogous procedure 
to the bottom mass. In particular, a threshold due to nib is added and one-loop SUSY 
corrections due to just chargino exchange are included after the initial step. 

At mz the initial values of ^ib(^z)susY an d m r(^z)susY are used to compute 
the corresponding Yukawa couplings h b and h T for a fixed value of tan (3. The Yukawa 
coupling h t of the top is computed at m t {m t ). These couplings are next run up to 
Mqut using two-loop RGEs in the DR scheme with the same initial mass threshold 
as for the gauge couplings, but are not assumed to unify. 

In the second step, starting from Mqut, for a given choice of tan /3, mi/2, mo, A Q 
and sgn(n), the RGEs for the masses and couplings are run down to the low scale. 
Here much attention is paid to extracting the parameters at a right energy scale 
Qq. The RGEs for the gaugino masses are run down until Q reaches their respective 
running values and analogously for the sfermion soft masses and the third-generation 
trilinear parameters A t , A b and A T . The top Yukawa RGE is run down to Qo = m\ ole 
while the RGEs for the other two Yukawa couplings of the third generation are run 
down to mz- 

Of particular importance is a correct treatment of the Higgs sector and the condi- 
tions for the EWSB. This is because of the spurious but nagging Q-scale dependence 
of the MS scheme 0]. Even including full one- loop corrections to the Higgs potential 



is not sufficient to significantly reduce the scale dependence when one minimizes the 
Higgs potential at the scale mz- Instead, it has been argued, for example in ||45| , that 
the scale dependence is significantly reduced by evaluating the (one-loop corrected) 
Higgs potential at Qj ~ y/ m t 1 m T 2 whh (mj 2 ) denoting the physical masses of 
the stops. This is because, at this scale, the role of the otherwise large log-terms 
~ log (m~/Q 2 ^J from the dominant stop-loops will be reduced. Here we follow this 
choice which is also adopted in ISASUGRA. (Actually, in ISASUGRA a very similar 
prescription is used with Qj= ^mj L mj R , where mj L (mj R ) denote the soft masses of 
the stops, but we do not think that this difference is of much importance). At this 
scale one evaluates the conditions for the EWSB which determine /i 2 as well as the 
bilinear soft mass parameter B\i 

, K + ^Q-K+^tan^ ! 
" tan 2 /3-l 2 mz ^ 
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2Bp = tan 2/3 [(m 



(m^ + S«)] + m|sin2/3 



(3.3) 



where m 2 H are the squares of the soft mass terms of the Higgs doublets B\^ u and 
E^u are their respective one-loop corrections. 

In the package ISASUGRA only the (dominant) third-generation (s)fermion one- 
loop corrections were included. We have added all the one-loop corrections following 
Ref. fr42], in particular the chargino and neutralino ones. It has been claimed in 
Ref. jE| that, while normally subdominant, these can contribute to increasing tua 
at large tan/3 and large m^j. 

The mass of the pseudoscalar will play a crucial role in computing Q x h 2 , es- 
pecially at very large tan/? ~ 50. This is so for three reasons: uia becomes now 



much smaller |49| than at smaller tan/3 due to the increased role of the bottom 
Yukawa coupling; because the A-resonance in XX ~~ * ff is dominant since the cou- 
pling Aff ~ tan /3 for down-type fermions; and because, in contrast to the heavy 
scalar H, this channel is not p-wave suppressed In ISASUGRA rriA is computed 
as 

m\ = (tan (3 + cot (3) (-Bfi + A^) (3.4) 

where A 2 A stands for the full one-loop corrections which can be significant |^2| . 

The neutralino relic density can be reliably computed both away from resonances 
and new final-state thresholds, where the usual expansion in powers of x = T/m x 
(where T is the temperature) works well, as well as in the vicinity of such special 



points. Exact analytic cross sections, which will soon become available [46], allow us 



to consistently and precisely determine Q x h 2 both near and away from such special 
points by solving the Boltzmann equation as in Ref. ||7|, (In the case of the 
expansion, analytic formulae for the thermally-averaged product of the annihilation 
cross section and the relative velocity are given in Refs. [T]] and are applicable far 
away from resonances and thresholds. However, caution is advised in combining them 
with a numerical integration of exact cross sections in the vicinity of the resonances. 
This is because one then usually neglects interference terms which can play some 
role. Furthermore, the range of m x around the pseudoscalar resonance where the 
expansion fails badly can be as large as several tens of GeV ||50|| .) 

We also include the co-annihilation with next-to-lightest SUSY particles (NLSPs) 
which in some cases is important. The co-annihilation with the lightest chargino 
and the next-to-lightest neutralino is treated without any approximation following 



Refs. [0, Hfl. In tne CMSSM, the lighter stau (n) is the LSP in the region m 1/2 > 



mo |T5[ and just above the boundary the neutralino LSP can efficiently co-annihilate 
with it and with the other sleptons p2[ . In this analysis we use the approximate 
expressions given in the second paper of Ref. |52| even though they do not include 
the effects of h T , of the T\ — t<2 mixing, in some channels of the mass of the r, etc 0], 
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which make them less reliable at large tan/5 > 20 j|. We have made an attempt 
at improving the available formulae by including in the propagators the widths of 
the gauge and Higgs bosons and the neutralinos, which otherwise become singular. 
Furthermore, ISASUGRA does not include one-loop corrections to the neutralino, 
chargino and slepton masses which can be of the order of a few per cent ||42| . This 
will have a comparable effect on the exact position of the boundary between the 
neutralino and the f\ LSP especially at large mx/2- For these reasons we will treat 
our numerical results in the co-annihilation region as somewhat less reliable but 
would not expect any major changes if the mentioned effects were included. 



4. Results 

We present our results in the plane (m^, m ) for several representative choices 
of tan/3 and other relevant parameters as specified below. First, in Figs. |T]-^| we 
show the experimental and cosmological bounds for respectively tan (3 = 10, 40 and 
50 and for A = 0, fx > and the central values of m t = mj? olc = 175 GeV and 
m b = m b (m b )fi = 4.25 GeV. In Figs. and | we can see many familiar features. 
At mi/2 ^> mo there is a (dark red) wedge where the fx is the LSP. On the other 
side, at mo mx/2 we find large (grey) regions where the EWSB is not achieved. 
Just below the region of no-EWSB the parameter /i 2 is small but positive which 
allows one to exclude a further (light red) band by imposing the LEP chargino mass 
bound. As one moves away from the wedge of no-EWSB, /i 2 increases rapidly. That 
implies that, just below the boundary of the no-EWSB region, the LSP neutralino 
is higgsino-like but, as one moves away from it, it very quickly becomes the usual 
nearly pure bino. This causes the relic abundance £l x h 2 to accordingly increase 
rapidly from very small values typical for higgsinos in the hundred GeV range, 
through the narrow strip (Ami/ 2 ~ 20 GeV) of the cosmo logically expected (green) 
range ( |1.1|) to much larger values, excluded (light orange) by (|1.2j ). In particular, in 
the whole region allowed by the chargino mass bound the LSP is mostly bino-like. 
We remind the reader that in the nearly pure bino case the neutralino mass is given 
by m x ~ 0.44m 1/2 - 2.8 sin 2(3 @. 

It is clear from Figs. [I] and that up until tan/3 ~ 40 the overall shape of 
the cosmologically expected ( |1.1| ) and excluded ( |L2p regions does not change much. 
Generally one finds a robust (green) region of expected Q x h 2 at ~ m in the 
range of a few hundred GeV fp~5]| . In addition, at ^> m , just above the wedge 
where the LSP is the fx, the co-annihilation of the neutralino LSP with f\ opens up a 
very long and very narrow strip which is allowed by the bounds ( |1.1|) - (|1.2|) . Finally, 
as mentioned above, at uiq ^> m^, very close to the region of no-EWSB, again one 
finds a very narrow range of Q x h 2 consistent with ( |1 . 1| ) . A closer examination would 
be required to establish to what extent this region of m 1( / 2 and tuq overlaps with that 



corresponding to 'focus points' advocated in Ref. [[20[ . 
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The region of no-EWSB is quite sensitive to the relative values of the top and bot- 
tom masses. Generally, at fixed tan/3, increasing (decreasing) the top mass relative 
to the bottom mass causes the region of no-EWSB to move up (down) considerably 
because of the diminishing (growing) effect of the bottom Yukawa coupling on the 



loop correction to the conditions of EWSB (p.2|)-(p73|). This effect can be explicitly 
seen in Fig. |4]a and b but it remains basically true also for much smaller tan (3. At 
fixed top and bottom masses, as tan (3 decreases, the region of no-EWSB moves to- 
wards somewhat larger values of tuq and smaller values of m\i<i but the overall effect 
is not very significant. 

The lightest Higgs mass also grows with mt as expected. In the Figures we plot 
the contours of mh = 113 GeV (which, within the uncertainties mentioned earlier, 
corresponds to the LEP lower limit) and of rrih = 115 GeV where an intriguing 
possibility of a Higgs signal has been reported ||21|| . Larger values of m^ are given 
by contours which are shifted along the m\ji axis with roughly equal spacings but 
diverge somewhat at larger mo- We remind the reader that the values of rrih that 
we plot have been obtained using ISASUGRA. When we use FeynHiggsFast p5 |, we 
obtain values lower by ~ 2 — 3 GeV, as can be seen from Table |l[ 

Other constraints behave with increasing tan/3 as expected. In particular, the 
(light brown) region excluded by BR(£> — > X s j) grows significantly because the 
dominant chargino-squark contribution to the branching ratio grows linearly with 
tan (3. On the other hand, the excluded region does not change much as we vary the 
assumed lower bound on BR(5 — > X a j). Taking the (slightly higher) most recent 
range (see Section 0) enlarges the excluded regions of the (m^, mo)-plane by some 
additional 50 GeV. 

As for the SUSY contribution to the anomalous magnetic moment of the muon, 
the lcr, and especially 2a, (yellow) range of allowed Aa^ USY becomes significantly 
enlarged towards larger m-1/2 and mo because the dominant sneutrino-chargino loop 
exchange contribution grows with tan/3 [54]. It is clear that the dependence of 
Aa^ USY on mi/2 and mo is relatively weak and therefore, in our opinion, no robust 
upper bound on superpartner masses can at present be placed. On the other hand, 
the low ranges of mi/2 and mo below the lines corresponding to the upper la (2a) 
limits on Aa^ USY = 59, 75 x 10 -10 , respectively, are probably firmly excluded but are 
not too big. 

The region of very large tan/5 ~ 40 — 50 deserves a separate discussion. First, 
as can be seen by comparing Figs. [I] and 0, we find relatively little variation in the 
overall shape of the constraints as tan/3 increases up to ~ 40. However, as tan/3 
grows, above ~ 45 a new important feature appears which is prominently reflected 
in the neutralino relic abundance, as can be seen from Figs. |3| and [|. The heavy 
Higgs bosons become 'light' enough j|9], [T^, |53], Q to "enter the (m-i/2, mo) -plane 
from the right" and cause a significant decrease in Q x h 2 , predominantly through a 
very wide A-resonance. Clearly, unlike for smaller values of tan (3 < 45, there is now 
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much variation in the position of the A-resonance and accordingly in the shape of 
the allowed and excluded regions of Q x h 2 as one varies m t and/or m&. Generally, 
as nib increases relative to m t , the position of the resonance moves to the left and 
the shapes of the regions where fi x h 2 is in the prefered range fll.lj) and is excluded 
by ( |1 . 2| ) change accordingly Qualitatively this behavior again clearly demonstrates 
the growing impact of hb at large tan/3. More quantitatively, one can see why this 
happens by examining Eq. ( |3.4|) where the one-loop correction A 2 A Jl6| to m\ from 
the sbottoms grows with tan/? as m\j cos 2 /? ~ hf. 

The white regions to the left of the green regions of 0.1 < Q x h 2 < 0.2, as well 
as around the green 'island' in Fig. (|d, correspond to £l x h < 0.1. The white areas 
between the favored green regions and the excluded light orange range (Q x h 2 > 0.3) 
as well as the white pocket inside the green region in Fig. ^d correspond to 0.2 < 
Q x h 2 < 0.3. One should also mention that, in contrast to much smaller tan/?, at 
large tan/3 ~ 50, Q x h 2 does not grow rapidly with increasing superpartner masses 
but instead changes very gradually and does not exceed unity by more than a factor 
of a few. This is again mostly caused by the very wide A-resonance because the A- 
width is large and actually growing with tan/? as well as mi/ 2 and m . As one cuts 
along the mi/2 axis at fixed m in Figs. ^| and £|, one can see that fl x h 2 first grows 
somewhat, then gradually decreases around the A- resonance and then increases again 
at large mi/2. 

The most recent measurements of Qu and ^cdm have now implied much more 
restrictive ranges 0.10 < Q x h 2 < 0.12 and 0.10 < Q x h 2 < 0.15, as discussed earlier. 
We do not show these ranges in our Figures but can easily summarize their effect. 
The range 0.10 < Q x h 2 < 0.15 corresponds roughly to the half of the allowed green 
strips (when cutting along them) on the side of the axes and the origin. In Fig. §1 this 
corresponds to the the outer halfs of the green 'island'. Requiring 0.10 < Q x h 2 < 0.12 
causes a further shrinking of the green regions in the same direction by roughly 
another factor of two. 

Clearly, as tan/? increases above ~ 45, the cosmologically expected and excluded 
regions are gradually shifted towards larger and m thus significantly relaxing 
the tight bounds characteristic of smaller tan/?. This is especially true when the top- 
to-bottom mass ratio is on the lower side. Overall, however, cosmological constraints 
on Q x h 2 now permit much larger superpartner masses than at smaller tan (5, and not 
only in the very narrow strips close to the regions of no-EWSB and/or ti-LSP. 

The case of very large tan /3 > 50 is obviously quite complicated. We are not even 
convinced that ISASUGRA and other currently available codes for generating SUSY 
mass spectra are fully reliable at such large values of tan /?, especially in the regime 
of mi/2,m > 1 TeV. First, two-loop corrections to the effective potential would 
need to be added to ensure a further reduction of the scale dependence. Second, now 
the Yukawa couplings of the bottom and the tau become comparable with that of 
the top and the contribution to the effective potential from the sbottom and stau 
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loops becomes of the same order as that from the stops. This introduces at least two 
more mass scales and the choice of the right scale for minimizing the Higgs potential 
becomes much more complicated. The variation in the shape and position of the 
region of no-EWSB and the heavy Higgs masses with m t and that we discussed 
above is, in our opinion, probably just a reflection of the above problems. 

One may conclude that, given the technical difficulty involved, the available codes 
for computing especially /i and the heavy Higgs masses are not yet fully applicable 
to such large values of tan/3, especially for m 1 /2,m > 1 TeV and we applaud the 



ongoing efforts to ameliorate the situation |56 

On the other hand, the case of tan/? ~ 50 is an intriguing one as it corre- 
sponds to the infrared quasi-fixed point solution of the top Yukawa coupling for 
the Yukawa RGEs when the approximate unification with and h T is further as- 
sumed |57[ [58], p9[ . Additionally, in S'O(10)-based models various texture models 
for mass matrices of leptons and quarks invariably require large tan/? ~ 50 |J0], ^Hfl . 
Collider phenomenology can also be distinctively different |62| with the usual mul- 
tilepton signatures for SUSY signals now being much reduced, although new signals 
may appear involving r-leptons and 6-quarks in the final state. 

Finally, in Table [3] we present several representative cases consistent with all 
experimental and cosmological bounds. All the mass parameters are given in GeV. 
The Table should be helpful in a more detailed comparison of our results with other 
groups. 

Our results show a sizeable difference at large tan f3 3> 10 when compared with 
some other recent analyses [p3 , |], First, in Ref. |k| a similar procedure for 



computing Higgs and SUSY mass spectra was used as in our modified version of 
ISASUGRA, which was based on Ref. 



421. In Ref. 



somewhat different quark 
input masses were used (m]? ole = 174 GeV and ml° le = 4.9 GeV) but, even for that 
choice, we find that our version of ISASUGRA still produced significantly lower 
values for rriA (by some 100 GeV) at large tan f3 ~ 50 than in Ref. |63 |, which we find 
somewhat puzzling. (This affects the position of the A-resonance and therefore also 
the relic density contours, even though in Ref. |)3| the latter was computed only in an 
approximate way using the package NEUTDRIVER.) As regards Refs. BJ, we find 
a sizeable discrepancy especially in the mass of the heavy Higgs bosons (especially 
the pseudoscalar) and therefore also in the position of the A-resonance. We believe 
that our computation of Q x h 2 , especially near the A-resonance, is somewhat more 
precise due to our exact treatment of the xx-annihilation. However, we are convinced 
that the main difference comes from the fact that in Refs. 0, || the effective Higgs 
potential was minimized at the scale Qo = mz while in ISASUGRA this is done 



at Q? 



as discussed in Sect. K This leads to large differences in the 



values of fi and B\i especially at larger ttiq and mi/2 which were also discussed in 
Ref. |64]]. We have checked numerically that choosing the minimization scale at mz 
typically gives significantly smaller mass of the pseudoscalar for a given choice of 
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Case 


RRN1 


RRN2 


RRN3 


RRN4 


RRN5 


RRN6 


RRN7 


tan p 


10 


10 


40 


45 


50 


50 


50 


mi/2 


300 


500 


250 


600 


700 


280 


1000 


m 


75 


1900 


1225 


350 


600 


1250 


1000 


pole 

ml 


175 


175 


175 


175 


175 


175 


175 


m b( m b)sM 


4.25 


4.25 


4.25 


4.25 


4.25 


4.25 


4.25 


A t 


-556.3 


-886.8 


-453.8 


-1028 


-1184 


-495.3 


-1644 


A 6 


-803.5 


-1245 


-573.5 


-1274 


-1412 


-579.2 


-1937 


A T 


-186.0 


-301.5 


-74.2 


-166.5 


-143.0 


-38.4 


-219.6 


m h (ISA) 


114.6 


118.0 


115.5 


120.0 


120.7 


115.8 


122.3 


m A (FHF) 


112.8 


115.1 


113.1 


118.3 


119.0 


113.5 


120.5 


m A 


443.6 


1920 


833.6 


601.9 


652.5 


537 


917.5 




397.3 


299.4 


166.7 


701.5 


785.2 


209.5 


1038 


m x 


117.5 


196.6 


87.8 


247.9 


291.7 


105.3 


424.6 


m xf 


215.5 


272.7 


129.6 


461.1 


542.4 


163.7 


784.6 


m~ 


706.3 


1224 


659.5 


1333 


1558 


727.9 


2176 


m e L 


221.0 


1919 


1231 


515.9 


757.2 


1258 


1193 


m e R 


140.9 


1905 


1227 


394.8 


653.2 


1253 


1064 


ra- 


132.2 


1889 


1035 


254.4 


429.6 


911.0 


734.1 


ni- 


224.4 


1911 


1140 


503.1 


695.6 


1104 


1075 


m u L 


630.0 


2085 


1294 


1219 


1477 


1341 


2118 


m u R 


608.9 


2075 


1293 


1180 


1434 


1337 


2049 


m d L 


634.9 


2087 


1296 


1221 


1479 


1343 


2120 




607.8 


2074 


1293 


1175 


1428 


1337 


2041 


468.3 


1398 


835.8 


948.0 


1139 


870.0 


1634 


m t 2 


651.9 


1780 


1021 


1137 


1335 


1020 


1869 


Tftr 

01 


580.3 


1780 


1003 


1049 


1255 


996.5 


1806 


mr 


609.6 


2059 


1135 


1123 


1331 


1093 


1877 


BR(B -> X s7 )/10- 4 


3.4 


3.6 


2.8 


2.6 


2.8 


2.6 


3.0 


a,/10- w 


20.5 


1.1 


12.0 


21.8 


15.0 


14.4 


7.0 


Q x h 2 


0.17 


0.20 


0.16 


0.14 


0.13 


0.17 


0.16 



Table 1: Masses and other variables for several representative choices of parameters. 
Quantities with mass dimension are given in GeV. The mass of the lightest scalar Higgs 
is computed using ISASUGRA (ISA) and the FeynHiggsFast (FHF) routines. 



other parameters. This has a strong effect on the location of its resonance in the 
(mi/2, m )-plane and therefore on the neutralino relic abundance. It also strongly 
affects the position of the region of no-EWSB. 
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5. Summary 



In summary, the interplay of the different experimental constraints and the cosmo- 
logically expected ( [Lip and allowed ( |1.2| ) regions shows a significant dependence on 
tan (3. Generally, the region excluded by BK(B — > X s ^y) grows with tan f3 and so 
does the 2a region favored by the recent measurement of the anomalous magnetic 
moment of the muon. The lightest Higgs mass limit rrih > 1 13 GeV gives a lower 
bound on m^, which however becomes weaker as ijiq and/or m t grow. Additionally, 
there still remains some theoretical uncertainty in computing rrih, of the order of 
2-3 GeV. As can be seen by comparing the rrih contours of 113 and 115 GeV, this 
affects the range of excluded by several tens of GeV at larger m t and even by 
nearly hundred GeV at smaller m t . 

At tan (3 < 40, the experimental bounds exclude large parts of the generic 'low- 
mass' region of both mi/ 2 and m (below a few hundred GeV [[TJJ) of the expected 
(green) range of the neutralino relic abundance (0.1 < Q x h 2 < 0.2) which itself is 
strongly squeezed by the conservative upper bound £l x h 2 < 0.3. What remains are 
the two narrow bands of expected VL x h 2 along the wedges of no-EWSB and fi-LSP 
of which the second is (slightly) more favored by the measurement of a M . 

The picture changes dramatically as tan/3 approaches 50. The cosmologically 
expected regions expand significantly both in size and towards larger m\/2 and mo. 
They also become further affected by the appearance of the resonance due to the 
pseudoscalar, whose mass decreases significantly at very large tan (3 > 45, and es- 
pecially for smaller values of the ratio of the top to bottom masses. Overall, it is 
clear that, in contrast to smaller values of tan/3, at such large tan/? cosmological 
constraints point towards significantly larger gluino and/or sfermion masses in the 
TeV range. 

NOTE ADDED: After our analysis was completed, we became aware of Ref. Jfj 
where the impact on the CMSSM of the recent highly restrictive determination of 
^cdm^ 2 — 0.13 ± 0.01 [|l2j was discussed, along with other constraints, except for 
b — > S7, for the choice m^ ole = 175 GeV and ^(m^g^ = 4.3 GeV p5| . We find 



a reasonable agreement in the overlapping results, in particular in the location of 
the resonance of the pseudoscalar A. Similarly, we find a reasonable agreement with 
Refs. [551 |Ew|| which were completed after our paper. 
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tanjS=10 m t -175 GeV, m b =4.25 GeV 




500 10O0 1500 2000 



m 1/2 (GeV) 



Figure 1: The plane {jni/2i m o) f° r tan/? = 10, Aq = 0, /i > and for mt = m-f ole = 
175 GeV and m& = mb(mb)^ = 4.25 GeV. The light red bands on the left are ex- 
cluded by chargino searches at LEP. In the grey wedge in the left-hand corner elec- 
troweak symmetry breaking conditions are not satisfied. The dark red region denoted 
*X NOT LSP' corresponds to the lighter stau being the LSP. The large light orange 
regions of ^ x /i 2 > 0.3 are excluded by cosmology while the narrow green bands corre- 
spond to the expected range 0.1 < £l x h 2 < 0.2. Also shown are the semi-oval contours of 
a^ USY = Aa^ Y /10 -10 favored by the anomalous magnetic moment of the muon measure- 
ment at 2a CL (aS USY = 11,75) and la CL (aS USY = 27,59). The 2a range is shown in 
dark yellow. The three lines in the figure correspond respectively to a^ USY = 11, 27 and 59, 
when moving towards the origin. The lines of the lightest Higgs scalar mass = 113 GeV 
and 115 GeV are denoted by short and long-dash lines, respectively. 



16 



tan/? = 40 m t =175GeV, m b -4.25 GeV 




500 1000 1500 2000 

m i/2 (GeV) 

Figure 2: The same as in Fig. [l] but for tan (3 = 40. In addition, the light brown region 
now appears which is excluded by the lower bound of BR(i? — > X s ^f) = (3.11 ±0.69) x 10~ 4 . 
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tan/S = 50 m t =175 GeV, m b =4.25 GeV 




m m (GeV) 

Figure 3: The same as in Fig. | but for tan/3 = 50. The green band of the expected 
range 0.1 < Q x h 2 < 0.2 has now changed considerably due to the appearance of a wide 
resonance XX ~ > A — * ff. The white areas closer (further away) from the axes correspond 
to n x h 2 < 0.1 (0.2 < n x h 2 < 0.3). 
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tan/S=50 



11^=170 GeV, m„=4.25 GeV 



tan/S = 50 



m t =180GeV, m b = 4.25 GeV 
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Figure 4: The same as in Fig. |3| but for different values of the top and bottom masses as 
indicated in respective windows. In window d) the white area around the green 'island' of 
the expected range 0.1 < r2 x /i 2 < 0.2 corresponds to £l x h? < 0.1. The white strips between 
the narrow green bands and the excluded (light orange) regions of S7 x /i 2 > 0.3 correspond 
to 0.2 < Q x h 2 < 0.3, and so does the white 'hole' inside the green 'island' in window d). 
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